Primitive turbulence: kinetics, Prandtl's 
mixing length, and von Karman's constant 

Helmut Z. Baumert 

Freie Universitat Berlin, Germany 

Dept. Mathematics and Computer Science 

Abstract. The paper presents a theory of shear-generated turbulence at 
asymptotically high Reynolds numbers. It is based on an ensemble of dipole 
vortex tubes taken as quasi-particles and realized in form of rings, hairpins 
or filament couples of potentially finite length. In a not necesserily planar 
cross sectional area through a vortex tangle, taken locally orthogonal through 
each individual tube, the dipoles are moving with the classical dipole velocity 
u — r ■ to. The vortex radius r is directly related with Prandtl's classical 
mixing length. The quasi-particles perform dipol chaos which reminds of 
molecular chaos in real gases. Collisions between quasi-particles lead either 
to particle annihilation (turbulent dissipation) or to particle scattering 
(turbulent diffusion). These ideas suffice to develop a closed theory of shear- 
generated turbulence without empirical parameters, with analogies to birth 
and death processes of macromolecules. It coincides almost perfectly with 
the well-known K-Q turbulence closure applied in many branches of science 
and technology. In the case of free homogeneous decay the TKE is shown 
to follow t" 1 . For an adiabatic condition at a solid wall the theory predicts 
a logarithmic mean-flow boundary layer with von Karman's constant as 
(2 7T)- 1 / 2 « 0.399. 

"Everything should be made as simple as possible, but not simpler. " 

- Albert Einstein 

Keywords: Turbulence, vortex dipoles, vortex gas, dipol chaos, quasi- 
particles, von Karman's constant, law of the wall 



Introduction 

" The diversity of problems in turbulence should not ob- 
scure the fact that the heart of the subject belongs to 
physics." (Falkowski and Sreenivasan 2006). Most clas- 
sical turbulence theories and models rest more or less 
on the Fridman-Keller (1924) series expansion of the 
Navier-Stokes equation, the first element of which is 
the Reynolds (1895) equation. Higher elements are the 
subject of turbulence closure models (e.g. Wilcox 2006). 
This approach attracted many researchers, but could 
not solve most elementary problems of turbulence. On 
the contrary, it led to more and more 'universal' higher- 
order equations. Although also the writer of this article 



made his first steps towards turbulence along this path, 
he became more and more sceptic that it is the right one 
to proceed. Therefore, with the exception of the mean- 
flow kinetic energy balance, the following text does not 
take much advantage from the Fridman-Keller expan- 
sion. 

Below a new path is testct. It is mathcmaticlly sim- 
ple, and for some readers even primitive so that the 
writer felt obliged to announce this already in the head- 
line. Nevertheless, a primitive model sometimes high- 
lights aspects which are overseen in more elaborate the- 
ories (Lorenz 1960). 

The following development explores physical analo- 
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gies between turbulence taken as an ensemble of 
vortex-dipol tubes and real gases of semi-stable macro- 
molecules. The underlaying non-linear statistical the- 
ory of macromolecular gas kinetics is not repeated here 
and can be found, for example, in textbooks on syner- 
getics (Haken 1978) . 

While in ideal gases molecules are point masses with 
zero cross section and infinite free path, in real gases the 
cross sections are finite so that molecule-molecule col- 
lisions limit the free path. A collision of two molecules 
(higher-order collisions arc neglected here for simplic- 
ity) can have one out of two possible results: 

(i) The molecules are spatially scattered, which cor- 
responds to molecular diffusion and mixing. 

(ii) The molecules vanish due to their semi-stability, 
which corresponds to particle annihilation or con- 
centration decay. 

The turbulent analogues of the above two processes are 
(i) turbulent diffusion and (ii) annihilation of vortices 
or dissipation of turbulent kinetic energy. The particle 
dynamics (molecules, vortices) can be described by the 
following equation, which is a special case of the Orego- 
nator (a certain class of diffusion-reaction systems, see 
e.g. Ch. 9 in Haken, 1978): 

dn d ( dn\ 2 

Here n is the particles' volume density, v is the diffu- 
sivity, and [3 is a constant. 

Equation describes the joint effect of two 
irreversible processes: particle diffusion through 
d/dx (vdn/dx), and irreversible decay through —fir?. 
It forms an initial- value problem for n = n{x,t). The 
final state is n(x, oo) = 0. Finite molecule numbers can 
only be kept when new particles are continuously added 
to the volume under study. 

The following theory uses discrete vortex dipoles in- 
stead of the a.m. semi-stable macromolecules. They 
are quasi-particles, i.e. locally excited states of an oth- 
erwise homogeneous fluid like phonons in solid-state 
physics. At very high Reynolds number Re we ex- 
pect that the fluid contains asymptotically many vor- 
tices and exhibits universal behavior in the sense of the 
Kolmogorov (1941) scaling laws. 

Vortices 

The classical closed vortex line (e.g. the closed center- 
line of an idealized smoke ring) represents an exact weak 



solution of the Euler equation, has infinitely thin diam- 
eter, infinitely high angular velocity, but finite circula- 
tion, c = 7r r 2 to < oo, where r and u> are effective radius 
and vorticity. The fluid outside the vortex line is invis- 
cid and irrotational. The classical vortex is an abstract 
frictionlcss object. 

The Batchelor couple (Lesieur 1997) is a vortex dipol, 
consisting of two anti-parallel vortex lines. It is also an 
idealized frictionless object. When isolated and far from 
boundaries, in 3D the trajectory of a Batchelor couple 
is a straight line. The flow field of one vortex moves the 
other vortex and vice versa. In practice such couples 
are stable over short to moderate propagation times. 
They conserve their kinetic energy, circulation (which is 
zero for this couple by definition) , their vortex radii and 
vorticities (which carry opposite signs). The Batchelor 
couple propagates into the same direction as the fluid 
moves between the two vortices forming the couple. In 
a sufficiently dense ensemble of Batchelor couples their 
trajectories are no longer straight lines due to mutual 
interactions. 

The counterpart of the Batchelor couple will be 
called below a von-Karman couple. It is a pair of par- 
allel vortex lines of non-zero total circulation. In an 
initial phase of their evolution they rotate around their 
common center of gravity, which remains in rest. Such 
a couple is known since long to be fundamentally un- 
stable. Its kinetic energy is dissipated into heat (e.g. 
Lamb 1932) . 

Vortex dipoles vs. vorticons 

A vorticon is a narrow relative of the Batchelor couple, 
with the following differences: the effective radius r 
and its kinetic energy u 2 /2 = r 2 w 2 /2 are finite. These 
conditions are well realized in practice, at best in 
quantum turbulence (cf. Vinen and Donclly 2007) . 

Vorticon generation = turbulence production: 

Due to the conservation principle for circulation in ideal 
fluids (Hclmholtz 1858), in a circulation- free volume 
vortices can be generated (and annihilated) only in form 
of anti-parallel vortex pairs with vanishing circulation. 
I.e. turbulence kinetic energy (TKE) generation by 
shear is generation of quasi-particles in the above sense. 

Vorticon annihilation = TKE dissipation: Vor- 
ticons can be annihilated in collisions if the two colli- 
sion partners are reorganized into von-Karman couples. 
Consider two initial vorticons, (+ ff — ) 1 and (— JJ. +) 2 , 
which reorganize during collision into the following von- 
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Karman couples: (+ 1 || + 2 ) and (— 1 || — 2 ), where the 
numbers 1 and 2 refer to the numbered particles. In 
both new 'particles' the non-zero molecular viscosity of 
a real fluid leads to kinetic energy dissipation and fi- 
nally to a merger of the like-signed vortices on a lower 
energy level (see also Klein & Majda 1993). This 'in- 
elastic' collision happens when the angle ip> between the 
propagation directions of the two colliding vorticons be- 
longs to the interval \ tt < ip < | tt (forward facing half 
sphere) . 

In the initial definition of our quasi-particles we have 
neglected viscosity. We re-introduce this important 
property by assuming that von-Karman couples are 
absolutely instable and instantaneously converted into 
heat. 

Vorticon scattering = turbulent diffusion: Vor- 
ticons are scattered in collisions between two vorticons 
if the two are reorganized into two new Batchelor cou- 
ples with new propagation angles. Consider again two 
initial vorticons, (+ ft -) 1 and (— JJ. +) 2 . After a col- 
lision they may form a new quasi-particle with a new 
propagation velocity vector: (+ 1 || — 2 ) and (— 1 || + 2 ). 
This 'elastic' collistion happens when the angle ip be- 
tween the propagation directions of the two colliding 
vorticons belongs to the backward facing half sphere. 

Vorticon tangle 

At high Re vorticons are assumed to form a dense three- 
dimensional tangle. As already found for 2D trajecto- 
ries (Aref 1983, Eckhardt and Aref 1988) , in 3D vor- 
ticons as elements of a tangle propagate along complex 
non-linear trajectories until collision with other vorti- 
cons. But we expect that the tangle is so dense that 
the effective trajectories are short and can be treated as 
straight lines. Frisch (1995) suggested that "The idea of 
conservative dynamics punctuated by dissipative events 
could be directly relevant in three dimensions." This is 
exactly what is elaborated here. Our form of 'conserva- 
tive dynamics' is chaotic dipol motion which conserves 
energy until particle annihilation. The trajectories are 
'punctuated' by dissipative collisions which lead to vor- 
ticon conversion into heat. 

Figure [l] exhibits a cross section through vorticons. 
The arrows denote the propagation direction. A plus 
sign indicates rotation against the clock. 

Historical notes 

Whereas the present approach to turbulence is based on 
ensembles of discrete quasi-particles, most turbulence 
closures of the past take the fluid as a continuum. Both 
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Figure 1. Two elementary vorticon- vorticon interac- 
tions. On the very left the reference vorticon is shown 
which is assumed to move to the right. The middle col- 
umn shows different orientations of the target vorticon 
which all lead to a dissipative (inelastic) collision, i.e. 
to annihilation of a quasi-prarticle. The right column 
shows potential target vorticons in elastic scatter po- 
sition corresponding to turbulent diffusion. A collision 
with these targets conserves energy. 

concepts are old. James Clerk Maxwell (1867) gave 
a detailed discussion of this dichotomy in connection 
with his kinetic theory of gases and traced it back to 
Democrit and Lucretius. Rene Descartes (1596 - 1650) 
spoke about 'tourbillons' forming the universe; based 
on the results of Herrmann von Helmholtz (1858) Lord 
Kelvin (1867a, b) coined the notion 'vortex atoms'; for 
a comprehensive overview see (Saffman 1992) . 

Most of later analytic efforts towards a kinetic theory 
of vortices were restricted to spatially two-dimensional 
cases wherein the Hamilton formalism is applicable 
(Onsager 1949) . Exceptions are vortex-based numeri- 
cal simulation receipes and rules. 

Marmanis (1998) was the first who proposed the vor- 
tex dipole as the fundamental quasi-particle of turbu- 
lence. He wrote: "The introduction of the vortex dipole 
chaos assumption permits one to derive a kinetic equa- 
tion for a 'gas' of vortex dipoles." He presented a me- 
thodical mixture of Onsager's theory with Bogoljubov's 
perturbation-kinetic approach with application to 2D 
inviscid turbulence. Dissipative 3D turbulence has not 



Baumcrt: Primitive turbulence 



4 



been treated by this author. 

Lions and Majda (2000) tried to develop a kinetic 
theory of 3D turbulence at high Re. They used system- 
atic asymptotic expansions of the Navicr-Stokcs equa- 
tion in the Fridman-Keller sense. On the one hand they 
tried to overcome limitations of the Klein-Majda theory 
(Klein and Majda 1993) , on the other hand they limited 
themselves to the quasi-parallel case where the Klein- 
Majda theory is constructively applicable through a rig- 
orous equilibrium statistical formalism. 

Chen et al. 2004 also aimed at the development of 
a physico-kinetic turbulence theory. They explicitely 
envisaged Boltzmann's approach and developed a com- 
plex formalism, but still without practically applicable 
predictions. 

Despite high efforts, these works did not yet lead to 
unique and constructive rules for the practical compu- 
tation of the eddy viscosity from mean flow data. Un- 
fortunately, with a few exceptions eddy viscosity is the 
most important turbulence parameter in computational 
aero-, hydro- and thermo-dynamics. 

Below turbulence balance equations are derived 
which describe the local-average dynamics of turbulence 
kinetic energy, K, mean vorticity, uj, and eventually 
eddy viscosity, v. These equations are only a specific 
format or a specific implementation of the primitive the- 
ory. In special cases they can be solved analytically. 

In more general cases one would try to solve them 
using finite-difference or finite-volume etc. methods. 
If the solution shall be approximated by the Monte- 
Carlo method then the writer recommends to start not 
from the continuum equations derived below, but bet- 
ter directly from the mechanical picture of the vorticon 
ensemble, preferably using the already existing eddy- 
collision methods where various vortex-filament primi- 
tives are already avaialable (Cottet and Koumoutsakos 
2000, Andeme 2008). 



tive vortex radii Tj and vorticities ujj. During the free 
flight time Tj the properties of the quasi-particle j are 
conserved. The particle's volume density is n/SV. The 
total TKE within SV is the sum of the kinetic energies 
of the individual vorticons: 
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Multiplication of (|2| |3| with SV gives 
K = SV-JCsv 

UJ = SV ■ 0J$v 
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Ksv, ^sv a re local volume densities of TKE and vortic- 
ity magnitude, respectively. They and K, uj are exten- 
sive variables by definition, i.e. they change when new 
particles with average properties k, uj are added to SV, 
as then K — ► (n + 1) k and u> — ► [n + 1) uj. 

k and uj are ensemble averages and as such intensive 
variables which do not change when new particles with 
average properties are added to the ensemble in SV: 



1 \ " 1 2 2 

jeSV 
1 , 
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Auxiliary quantities 



We introduce useful auxiliary quantities: the ordinary 
vorticity frequency, = 1/T and the often used con- 
stant k: 



Primary variables in primitive 
turbulence 



UJ 
K 



= 2nn = VL/h 



(2 7T) 



-1/2 



0.399 



(8) 
(9) 



TKE and vorticity 

For the construction of the eddy viscosity, which is pro- 
portional to TKE times unit time, and of the TKE 
dissipation rate, which is TKE per unit time, two con- 
stituing variables are necessary: TKE and a measure of 
time. Here the below-defined quantities K and uj are 
chosen. They are introduced now. 

Consider a small volume element SV populated by 
an ensemble of j = 1 . . . n dipoles with individual effec- 



Further below, k appears to be von Karman's constant. 
While \ujj\ = 2ir/Tj is an angular frequency, f2 is an 
ordinary frequency. 

Derived variables 

Effective radius R. R 2 is defined here as the weighted 
ensemble mean: 



R 2 = 



(E 
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It depends (inversely) on n and is thus an extensive 
variable. In freely decaying homogeneous turbulence 
(Section ) only the particle number decreases while the 
ensemble mean properties k and Q remain constant in 
time such that 



n 7r R 



const. 



(11) 



This reminds of a 2D version of the equation of state 
of a real gas telling us that the (planarized) locally or- 
thogonal cross sectional area through a dense vorticon 
tangle is compactly filled with dipoles. 
Dissipation rate. The rate of TKE dissipation, e, has 
the units of TKE per time [m 2 s -3 ]. Up to a factor £ it 
is governed by the product K ■ Q, : 



„ KQ „ kio n 
C = C — n 2 

7T 7T 



(12) 



e is an extensive variable. Below we show that £ = 
1. In freely decaying homogeneous turbulence (12) is 
identical with the quadratic term in such that in 
this case 

(3 = (koj/ir = const. (13) 

Eddy viscosity. Ludwig Prandtl (1925) noticed two 
simultaneous properties of his mixing length, namely 
that it (a) "...may be considered as the diameter of 
the masses of fluid moving as a whole in each individ- 
ual case; or again, as the distance traversed by a mass 
of this type before it becomes blended in with neighbor- 
ing masses . . ."; and also that it is (b) ". . .somewhat 
similar, as regards effect, to the mean free path in the 
kinetic theory of gases ..." (Bradshaw 1974) . 

We interpret Prandtl 's mxing length in terms of 'our' 



R in ( 10 1 because it is proportional to (a) the local mean 



radius of vorticons and, simultaneously, (b) to the free 
path. 

To explore this idea further we follow Albert Ein- 
stein's (1905) theory of molecular diffusion in fluids and 
use his expression for the coefficient of diffusion now for 
our coefficient of turbulent diffusion, or eddy viscosity 
v, in terms of a mean free path, A, and a mean free 
flight time, r: 



= A 2 /2r 



(14) 



For a vortex dipole these parameters can directly be 
given as follows wherein u is the mean (linear) advection 
velocity of a vortex dipole: 



(15) 
(16) 




For a quasi-steady turbulent vorticon tangle at Re = 
oo and far from solid boundaries the mean free path 
will clearly be very limited but it cannot vanish. The 
packing density should be dense but small enough that 
during a flight a particle can be replaced by a newly 
generated. This implies that (Baumert 2005b) 



A = 2R. 



(17) 



The free area of a free vorticon far from boundaries is 
(Fig. [2) 



A fr 



{ARf 



(18) 



We finally get the eddy viscosity as an intensive variable 
as follows: 



= K/nn = 2K/u = u) R 2 



(19) 



This equation is called the Prandtl-Kolmogorov rela- 
tion. 

Free homogeneous decay of turbulence 

For constant properties k, ui of our quasi-particles, fol- 
lowing ([6] [7]) the variables K and f2 are proportional 
to the vorticon density n and have therefore to satisfy 
the same diffusion-reaction equation (JTJ like n, because 
these properties are 'solidly mounted' at the dipoles. In 
the case of decay we have thus the following: 



K _ LU 
k Lit 



In analogy to ([lj it holds that 
dK 



dt 

dn 

~dt 



V0 K K z = 0, 

+ (3 n n 2 = 0. 



The reader easily sees that at large t 

K(t) = (fat)- 1 , 
fi(t) = (Pnty 1 , 



and 



v(t) 



0c. 



K(t) 



const. 



A" 



(20) 

(21) 
(22) 



(23) 
(24) 



(25) 



Equation ( 23 1 coincides with the results of a fairly gen- 



eral similarity analyses of the Navicr-Stokes equations 
by Oberlack (2002) and with the experimental results 
of Dickey and Mellor (1980). 
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Figure 2. Local cross section through a dipol tangle 
far from boundaries with maximum possible density at 
Re — ► oo. The dark grey vorticon in the left half square 
labelles the instantaneous situation, the 'empty vorti- 
con' in the right half square indicates the vorticon after 
its free flight from the left to the right into its collision 
position, where it is either annihilated or scattered. Far 
from boundaries a vorticon obviously 'needs' a cross sec- 
tional area of Af re e = (4-R) 2 )- 



Forced inhomogeneous turbulence 

In contrast to Section now the assumptions of free de- 
cay and spatial homogeneity are given up. This implies 
that we need to add source (Fk and Fq) and diffusion 
terms to the right-hand sides of equations (21 22 1: 



dK 

~dt 

on 
~dt 



(3 K K 2 
- (3q n 2 



= F, 



K 



d f dK 

dx 



F a 



8 



dfl 



dx dx 



(26) 
(27) 



The factors (3k and (3q are not longer expected to be 
constants. 

TKE production by mean-flow shear: Fk 

The specification of the TKE production term Fk for 
a shear flow can be obtained from textbooks (Wilcox 
2006, Schlichting and Gersten 2000), 



F K = vS 2 = 



K 



S 2 



(28) 



and Ui is the i th component of the mean flow velocity 
vector (U 1 ,U 2 ,U 3 ) T . 

Vorticity generation by mean-flow shear: Fq 



The above formula ( 28 1 for dK/ dt — Fk is classical hy- 
drodynamic knowledge. However, the specification of 
the vorticity generation term Fq is less trivial. Here one 
has to consider a fundamental argument of Henry Ten- 
nekes. It has been carefully discussed in the past, see 
Tcnnekes 1989; Baumert and Peters 2000, 2004; Kan- 
tha 2004, 2005; Kantha et al. 2005; Kantha and Clayson 
2007. We cast Tennekes' argument into mathematical 
form. He hypothesized that, on dimensional grounds, 
the length scale (here: L or R) cannot depend on the 
ambient shear for a neutrally stratified homogeneous 
shear flow. Since shear production Fk involves shear, 
Fq needs to be constructed such that the role of shear 
vanishes in the evolution equation for the length scale. 



The latter is generally derived from equation (10 1 as 



1 dR 2 

W ~dT 



1 dK 

K ~dt 



1 duj 
u> dt 



We insert dK/dt from (|28f in Q and find 

1 duj 



1 dR 2 

W ~dT 



S 2 ^ 



u> dt 



(30) 



(31) 



In terms of the present theory, Tennekes' hypothesis 
means dR 2 /dt = 0, which implies that 



= 2- 



1 dto 
~dt 



It is left to the reader to see that ( 32 ) gives 



dn 

~dt 



2tt 



(32) 



(33) 



As far as ( 10 1 is a pure ensemble-mean definition which 



considers neither diffusion nor annihilation of particles, 



the time derivative in ( 33 1 needs to be understood as 



the pure generation term in (|27|, Fq: 

Fq 



52 
2tt 



(34) 



The multipliers (3k and (3f, 



Let us summarize (26 27) and (28 34 1 as follows 
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where v is known as function of K and f2 through (19 1. Logarithmic law of the wall 



It remains to determine the still unknown multipliers 
(3k and (3n- This can be accomplished as follows. 
We note that the term (3k K 2 in ( 35 1 is identical with 



We insert (46 1 into the stationary form of (W2| so that 



we have to solve the following equation for S — S(z): 



the dissipation rate ( 12 1 of TKE 



2K d (\ dS 



= (3k K z 



(37) 



From the second equation in ( 37 1 we conclude with ( 19 1 
that 



2£- 
The solution is 

S(z) = 



1 dz \ S dz 



S 2 



dU _ V2iV(2C-l) 



(3k = 



c 



dz 



(47) 



(48) 



while from ( 25 ) it follows that 

(3q = v it (3 K ■ 



We insert ( 38 1 into ( 39 1 and see that 
(3n = -. 



(38) 



(39) 



(40) 



Integration of (|48j) gives the logarithmic law of the wall. 
In boundary layer theory the bottom shear stress is 
mostly defined in terms of the squared friction veloc- 
ity u 2 f , 



dU 

dz 



S. 



and with ( 46 1 it follows that 



Inserting (|38| |40| in (|35j |36| we get, after some algebra This allows to rewrite (|48l as follows 



using (19 1, an almost complete system wherein only £ 



is still to be specified: 
dK 



S(z) = 



dU 

dz 



K Z 



dK 
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u(S 2 



dx \ dx 



7T \2 



c^ 2 



(41) 
(42) 



with ii defined through 



«^/(2C-i)v/c. 



Turbulent boundary layer 

Consider a stationary boundary layer close to a plane 
solid wall at z — where z is the only coordinate of 
interest here. It points orthogonal from the wall into 
the fluid. The mean flow is parallel to the wall, i.e. 
Ui = U > and U 2 = U 3 = 0. Consequently with ([29 1 
we have to write 



Integration of (51 1 provides us with 

u f 



U(z) 



In 



(49) 
(50) 

(51) 

(52) 

(53) 



Mixing length L 

Consider the definition of the effective vorticon radius 



through ( 10 1. We solve this equation for K and express 



S = \dU/dz\ 



(43) 



the TKE in terms of R and ft as follows: 
K = 2n 2 R 2 tt 2 . 



(54) 



The diffusive TKE flux into the viscous sublayer at z = 
zq — * has to vanish, 



or 



{vdK/dz) z=z<) = 0, 
K(z) = Kq = const. 



such that in the stationary case (41 1 gives 



n = s/y/c. 



(44) 
(45) 

(46) 



In practical analyses of boundary layer turbulence, 
Prandtl's (1925) concept of mixing length is still often 
applied. It relates the eddy viscosity to the shear. Fol- 
lowing Hinze (1959, p. 279, eq. 5-2) in present notation, 
Prandtl defined his mixing length L as follows: 



v = L 



dU 



dz 



= L 2 S. 



(55) 



Due to our eddy viscosity formula (19 1 relation (55 1 
gives 



L 2 = 



K 
ttCIS ' 



(56) 
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so that in the neighborhood of a solid wall we get with 



(46 1 the following result: 



K = tt L 2 ft 2 y Q . 



Equating (54 1 with (57) gives 

R 



L 



C" 



1/4 



(57) 



(58) 



The physical meaning of L is easily understood by a 
glance at Fig. [3j If we may set ( = 1 then (58 1 gives 
together with the definition of n in ^ the following, 



L 2 = 2{ttR 2 ) , 



(59) 



which is the area of a vorticon, fully compressed into the 
form of a square of side length L. This is the asymp- 
totically maximum deformation which justifies to set 

C = i. 

While Prandtl's mixing-length concept was appli- 
cable only in the vicinity of solid boundaries so that 
it attracted respectful criticism (e.g. Kundu 1990, p. 
457; Wilcox 2006), our concept is a generalization of 
Prandtl's concept and works also far from boundaries, 
even in the free stream of stratified fluids where L 
may approach the Thorpe scale and/or the Ozmidov 
scale, depending on the conditions (Baumert and Pe- 
ters 2004). 

This section shall be closed with a summary of for- 
mulae relevant at solid boundaries, where z the distance 
from the wall: 



v 
L 

z 



u f L, 

K Z , 

L/k = R/k 2 . 



(60) 
(61) 
(62) 



Boundary compression 

Far from boundaries the dipole chaos is isotropic and a 
vorticon's cross section exhibits circles (Fig. [T]). Closer 
to a solid wall symmetry is increasingly broken and we 
see ellipses instead. Right at the boundary the ellipses 
are even deformed to an equivalent square (the asymp- 
totic case) with side length L = (2 tt i? 2 ) 1 / 2 , Fig. [3] This 
is the minimum area a vorticon can asymptotically at- 
tain: 

A bound = L 2 = 2irR 2 . (63) 

Notice that the fluid within the vorticon is free of fric- 
tion so that circle, ellipse and square are energetically 
identical. 




///\\ < 



///\\ 



Figure 3. Cross section through a vorticon sheet at 
a solid boundary. In this Figure we have r = R, i.e. 
the sketched quasi-particle is to be understood as an 
ensemble average. 



To estimate the degree of boundary compression of 



vorticons quantitatively we compare area ( 18 ) 'occu 



pied' by a free vorticon far from boundaries with the 



corresponding area ( 63 1 for a compressed vorticon at 
the very boundary: 



•A free 
•Abound 



(4R) 2 
2L 2 



16 R 2 
A-kR? 



4 

7T 



1.27. 



(64) 




/// \\ Solid boundary ///\\ 



Figure 4. Squeezed vorticons like in the dark-grey 
square of Fig. [3] The dotted lines symbolize the cir- 
culation pattern within the squeezed quasi-particles. 
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Discussion 

Kinetic theory and the k-u closure of Wilcox 

The specification £ = 1 closes the primitive turbulence 
theory. In compact form it reads as follows: 




These equations are structurally fully identical with the 
k-u) closure model discussed by Wilcox (2006) , which 
is in many aspects the best out of more than 10 other 
semi-empirical two-equation models. 'His' k and our K 
are equivalent. The eddy viscosity definitions of Wilcox 
and the present paper agree effectively up to a factor of 
0(1). 

However, 'his' u> (he calls it somewhat uncertain a 
'specific dissipation rate') differs significantly from our 

n. 

cu = ft/V/?*- 3.33 • (68) 

where (i* = 0.09 is one out of 6 empirical parameters of 
the Wilcox model. It reproduces the logarithmic bound- 
ary layer and implies k = 0.41. 

A major quantitative difference between the two 
models lays in the free homogenous decay. While both 
exhibit K <~ t~ m , we have m = 1, which agrees with 
Dickey and Mellor's (1980) high-i?e laboratory experi- 
ments and with Oberlack's (2002) theoretical result for 
the free decay in the Navier-Stokes equation. Wilcox' 
model is tuned to m = 6/5 = 1.2, which agrees with a 
number of other laboratory results and with Oberlack's 
(2002) result for the free decay in the Euler equation. 

Today it is not clear why many decay experiments 
lead to m > 1 . Possibly it is a matter of initial condi- 
tions, see e.g. Hurst and Vassilicos (2007): at high Re 
viscosity is comparatively small so that its regularizing 
effect on the initial spectrum towards a fully self-similar 
decay spectrum will take more time than at lower Re. 
In some cases this time may exceed the lifetime of tur- 
bulence. 

Von Karman's number 

The primitive theory makes a quantitative prediction of 
von Karman's constant. This allows a comparison with 
observations, experiments and competing other theo- 
ries. While the latter are rare, there is a substantial 



literature discussing the precise value of k and its error 
bounds. The controversy whether the logarithmic or 
the power law is applicable to boundary layers (Baren- 
blatt) shall not be picked up here. The controversy on 
the universality of n (Chauhan et al. 2005) is possibly 
only a matter of proper definition. If we restrict the def- 
inition of k to the case of favorable pressure gradients 
only (cf. Chauhan et al.) then we would get a unique 
solution which is supported also by the Chauhan et al. 
data: n = 0.40. 

A theoretical effort by Telford (1982) shall be men- 
tioned, who derived a formula which connects k with 
two semi-empirical parameters, the entrainmcnt and 
the decay constants. He predicts n = 0.37. Further 
theoretical work has been done by Yakhot and Orszag 
(1986) who found k = 0.372 based on renormalization- 
group approximations. Bergmann (1998) claimes that 
k = 1/e w 0.368, based on a new physical interpreta- 
tion resting on semi-empirical considerations. 

From an analysis of observational data Hogstrom 
(1985) concludes 0.40 ± 0.01, which is still accepted 
today among most practitioners. But Telford and 
Businger (1986) draw his analysis into question and 
underline the enormous methodical problems in con- 
nection with the evaluation of the experimental data. 

Another group of authors supports the lower values, 
e.g. Jinyin et al. (2002). They report many different k 
values obtained in the past from different experiments 
and simulations, e.g. back to Nikuradse, Laufer, Lawn, 
Perry and Abell, Long et al., Millikan, and Zagarola 
and Smits (loc. cit. Jinyin et al.). They summarize that 
the lower and upper bounds of k are 0.36 and 0.45, re- 
spectively, whereby the range of experimental Reynolds 
numbers involved is also quite broad. In this sense Za- 
noun et al. (2003) support k = 0.37 while Nagib et 
al. (2004) prefer 0.38, underlining the influence of non- 
stationary conditions. Andreas et al. (2006) report a 
somewhat higher 0.387 ±0.003 based on experimental 
material from the atmospheric boundary layer. How- 
ever, in the Princeton superpipe a value of 0.421 has 
been found for smooth flow at Re w 3- 10 7 (Smits 2007), 
too. 

Summarizing the hot debates of the past, a review 
by Jimenez and Moser (2007) on wall turbulence states 
boldly: "The Karman constant k rs 0.4 is approxi- 
mately universal." In practice this can hardly be dis- 
tinguished from our primitive-theoretical value 0.399 for 
Re — > oo. 

Nevertheless, doubts remain and the discussion by 
Talamelli et al. (2009) is justified. They underline the 
limitations of the various existing superpipe facilities 
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and propose a new one that will allow good spatial res- 
olution even at very high Reynolds number. The new 
facility is the international effort CICLoPE and under 
construction in North Italy. 

Kinetic theory and spectra of turbulence 

The above theory makes a precise prediction of von 
Karman's constant k. Would it be possible to use the 
present kinetic picture of turbulence also to determine 
asymptotic values for the Kolmogorov constants Ck 
and C w in the spectral energy distributions over fre- 
quency and wave number? With great probability the 
answer is no. Such an approach would demand that the 
lower bounds X u and ui u in the energy integrals could 
be specified: 

/>oo />oo 

K = C K e 2/3 / fc~ 5/3 dk = C u e oj- 2 duj . 

J X u J uj u 

(69) 

These bounds are not simply connected with the 
ensemble-mean values of vorticity and length scale. 

But it cannot be excluded that a dissipation theory 
of quasi-particles can be developed in analogy to a ra- 
dioactive decay chain, where large quasi-particles are 
broken by a collision into smaller ones etc. down to the 
(infinitely small) Kolmogorov scale where they are an- 
nihilated. This could be imagined as a cup with a tiny 
whole. It would keep all quasi-particles of size > in 
and let only those out which have zero size. Possibly 
this approach could find some attention. 

Conclusions 

The primitive picture of turbulence makes a minimum 
use of the Navier-Stokes equation. Vortex properties are 
derived from the Eulcr equation and dissipaton is intro- 
duced through an analogy to a birth and death process. 
Methodically our position is a neighbor of Maxwell's 
kinetic theory of real gases: 

• Real gas: Free molecule flights and their colli- 
sions are strictly governed by inertia and short- 
range forces of a non-stationary Coulomb field 
with moving point sources. Maxwell neglected the 
details of these interactions. He replaced them 
with geometrical and probability considerations 
and a statistical equilibrium assumption. He got 
the probability distribution of the molecular en- 
ergies whereby advective, rotational and potential 
energy of the molecules could be treated as sta- 
tistically independent. In an isolated real gas to- 
tal energy and molecule number are conservation 
quantities. 



• Turbulence: Here neither energy nor particle 
number are conserved, but the flux of mechan- 
ical energy into a volume clement is in quasi- 
equilibrium with the heat flux across the bound- 
aries. In vortex dipoles translation and rota- 
tion are coupled and statistically not indepen- 
dent. Nevertheless one may apply Maxwell's cen- 
tral idea here as well. Although the movements 
of our quasi-particles are strictly governed by 
well-known rules (Navicr and Stokes), we take 
them as unpredictable and replace them with 
geometry and probability arguments. Following 
Eward Lorenz (1960), such an approach stands 
today on much harder theoretical grounds than 
in Maxwell's times. 

We used only the following assumptions: 

1 . A turbulent fluid at Re — oo is a volume saturated 
with chaotically moving vortex-dipol tubes which 
are treated as quasi-particles. 

2. A new-born quasi-particle is equipped with en- 
ergy lost in from the mean flow formulated with 
the eddy- viscosity hypothesis, —(u'w 1 ) — v S. Its 
vorticity is governed by Tennekes' hypothesis. 

3. Advection of dipol tubes as quasi-particles follows 
simplest classical laws. 

4. If quasi-particles collide they are either scattered 
(50 %) or annihilated/converted to heat (50 %). 

5. If scattered they perform a Brownian motion with 
Einstein's diffusivity in Prandtl-Kolmogorov for- 
mulation. 
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